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ABSTRACT 

Aims. We study the timing and spectral properties of the intermediate polar MU Camelopardalis (1RXS J0625 18.2+733433) to determine the 
accretion modes and the accretion geometry from multi-wavelength, multi-epoch observational data. 

Methods. Light curves in different observed energy ranges (optical, UV, X-ray) are extracted. The timescales of variability in these light curves 
are determined using Analysis of Variance. Phase-resolved X-ray spectra are created with respect to the most prominent detected periodicities 
and each fitted with an identical model, to quantify the differences in the fitted components. 

Results. The published tentative value for the spin period is unambiguously identified with the rotation period of the white dwarf. We detect 
a distinct soft X-ray component that can be reproduced well by a black body. The analysis of data obtained at different epochs demonstrates 
that the system is changing its accretion geometry from disk-dominated to a combination of disk- plus stream-dominated, accompanied with a 
significant change in brightness at optical wavelengths. 
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1. Introduction 

Intermediate Polars (IPs) are close binaries that consist of a 
magnetic white dwarf and a Roche lobe-filling late-type main- 
sequence star. Mass loss from the secondary star occurs via the 
inner Lagrangian point, and the mass is accreted by the white 
dwarf. In contrast to the true Polars, the white-dwarf rotation is 
unsynchronised with the orbital motion but is more rapid due 
to the angular momentum gained from the accreted matter. 

In the case of a high accretion rate or a weak magnetic field, 
the magnetosphere will be smaller than the circularisation ra- 
dius and an accretion disk will form. The inner disk is truncated 
by the magnetic field of the white dwarf, which channels the 
matter via accretion curtains to the magnetic poles. Accretion 
arcs are formed along the footpoints of accreting field lines. 
Along these arcs, the accretion rate is assumed to vary. The 
heated plasma in the accretion zone is a prominent site of hard 
X-ray radiation, which due to the lighthouse effect is modulated 
on the spin period, co. Because of irradiation, the inner-disk rim 
may emit radiation that is modulated on this period. 

For stronger magnetic fields, some or all of the material 
from the infalling stream might be directly captured by the 
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field without going through an intermediary accretion disk. In 
the case of this stream-fed accretion, material transferred to the 
white dwarf retains the orbital motion of the secondary star, 
producing periodic variability of the beat frequency a>-£l, or 
other sidebands with the orbital frequency Q. 

Most IPs accrete predominantly by means of an accretion 
disk, V2400 Oph being a rare example of a sys tem with pure 
stream-fed accretion dHellier & Beardmorell2002l) . A number 
of IPs show evidence for accretion by means of both mech- 
anisms simultaneously: via accretion curtains beginning at 
the inner disk rim, and via so-called disk-overflow when a 
part of the ballistic stream bypasses the accretion disk out- 
side the orbital plane. In some systems repeated changes be- 
tween the differen t accretio n states were obs erved (TXCol, 
Norton etaf]ll997t FO Aqr. lEvans et allEooH) . Principal evi- 



dence for stream-overflow in IPs has been based on indirect 
means, i.e. the detection of certain sideband frequencies in 
their power spectra, which is not always unambiguous. For in- 
stance, sideband frequencies may also be produced by spin- 
pulsed light reprocessed in orbitally-fixed structures (a problem 
mainly in the optical and UV); alternatively, amplitude modu- 
lation of the orbital period (e. g. due to X-r ay absorption) may 
introduce spurious side-bands ( Warneil l986). 

An article by [Patterson ( 1994b introduces IPs in general, 
whereas iHellierl (120021) discusses the mechanisms of stream- 
fed accretion. 
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Fig. 1. The light curves of MU Cam, obtained simultaneously on March 31, 2005, by XMM (PN and OM) and with the AIP 
70cm-telescope (R). The data shown here are binned in time, 60s for the first, second, and forth panel (counted from top to 
bottom), 120s for the lower two. The R-band data are unbinned. 



The X-ra y source MU Camelop ardalis (the new name 
according to iKazarovets et al. d2006l) for the source for- 
merly known as 1RXSJ0625 18.2+733433) was id e ntified 
as an intermedi a te pol ar by lAraujo-Betancor et al.l (120031) 
and Istaude et"aL (2003). Optical photometry showed strong 
variability on mainly two timescales, which were tentatively 
identified with the orbital period of the system and the spin 
period of the white dwarf, P orb = 16987(23)s = 0.19661(27)d 
and P snin = 1187.2 46(4) s = 0.01374127(5)d, respectively 



Table 1. The observations analyzed in this work 



(Sta ude et al.l 120031) . It remained an open question whether 
the measured spin period provided the white-dwarf rotation, 
if it was twice the true rotation value, or if it was a side-band 
period with the orbital variation. 



Date (dd/mm/yyyy) 


Telescope/Instrument 


Duration 


28/03/2005 


AIP 70cm / R 


5.0h 


30/03/2005 


AIP 70cm / R 


8.6h 


20/04/2005 


AIP 70cm / R 


4.7 h 


31/05/2005 


AIP 70cm / R 


7.6h 


31/05/2005 


XMM / PN+MOS+OM 


12.5 h 


06/04/2006 
08/05/2006 


XMM / PN+MOS+OM 
AIP 70cm / R 


10.6 h 
4.4 h 



MU Cam is one of a handful of IPs that show a promi- 
nent soft X-ray com ponent additional to the hard component 
(IStaude et al . 2003). An XMM-Newton observation was pro- 
posed to identify undoubtedly the spin period of the white, and 
to determine the origin of the soft X-ray radiation. 

We report on X-ray and UV observations of MU Cam with 
XMM-Newton on both March 3 1 , 2005, and April 6, 2006, and 
coordinated ground-based optical photometry. 



2. Observations 

A short summary of all observations analyzed in this work is 
provided in Table Q] 

2.1. Data from XMM-Newton 

The X-ray counting EPIC instruments were operated in small- 
window mode with the thin filter. The EPIC-RGS devices were 
used in the standard configuration. The optical monitor (OM) 
was kept in fast mode with filter UVM2. 
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Fig. 2. The light curves of MU Cam, obtained on April 6, 2006, by XMM (PN, MOS1, MOS2, and OM). The data shown here 
are binned to 120s. 



On March 31, 2005 (revolution 972), an observation with 
XMM-Newton (observation ID: 0207160101), of a scheduled 
duration of 49ksec, was performed. Due to the high back- 
ground, both EPIC-MOS detectors were switched off after 
13.6ks, while the EPIC-PN observed for 30ksec. The OM was 
in operation for 45 ksec. Ground-based observations were per- 
formed synchronously in the R-band. The X-ray observations 
corresponded to 1.75 binary orbits and 25 white-dwarf spin cy- 
cles, in which the EPIC-PN collected ~3 7,000 source counts . 



A preliminary analysis was presented bv lStaude et al. I d2006h 



The observed mean source count rates were 1 .54 s 1 in the PN, 
0.38 s" 1 and 0.43 s 1 , in the MOS1 and MOS2, respectively; 
2.3 s- 1 in the OM/UVM2; and 0.031 s" 1 and 0.040 s- 1 , in the 
first order of RGS 1 and RGS2, respectively. 

In revolution 1158 (on April 6, 2006, observation ID: 
0306550101), XMM-Newton observed the target again for a 
scheduled ~39ksec. All instruments were used for the com- 
plete observation duration. Due to a telemetry gap, half an 
hour of data is missing for all instruments. Additionally, two 
blocks of 4400 s exposure in the OM were also lost. During 
this observation, the mean source count rates were measured 
to be 1.17 s- 1 in the PN; 0.30 s" 1 and 0.34 s- 1 in MOS1 and 
MOS2, respectively, 4.2 s-' in the OM/UVM2; and 0.022s- 1 
and 0.029 s 1 in the first order of RGS1 and RGS2, respec- 
tively. 



The data reduction was completed using XMM-SAS ver- 
sion 6.5; in particular the tasks epchain, emchain, and rgsproc 
were adopted to process the X-ray data. Source and background 
photons were extracted in adjoining circles for the PN and for 
the MOS data from 2006, while the background in the MOS 
observation of 2005 was extracted from an annulus around the 
source circle. Events with detection patterns of up to quadru- 
ples were selected. 

Background-subtracted light curves in different X-ray en- 
ergy bands - created with evselect - of the PN data are shown 
in Figs. Q] and [2] with the OM UV-data and ground-based R- 
band data acquired contemporaneously. The energy bands were 
selected after inspection of the shape of the averaged spec- 
tra (Fig. |4]i, which showed a soft component below 0.7 keV, 
and a minimum in the harder energy range between 2 and 
3 keV. We therefore define the following PN energy bands: 
SX: 0.3 - 0.7keV,MX: 0.7-2.0keV, HX: 2.0- 12.0keV, MHX: 
0.7 - 12.0keV, and TX: 0.3 - 12.0keV. 

EPIC source spectra were created with the task especget. 
RMF files for the data from the RGSs were created with rgsrm- 
fgen. In further analysis, all spectra were grouped with grppha 
to contain at least 20 counts per bin. 

The OM light curves were extracted using the task om- 
fchain. 
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2.2. Optical photometry 

Photometric observations in the optical wavelength range were 
performed with the 70cm telescope of the AIP in Potsdam. 
Differential photometry for this data was completed using 
DoPhot. The apparent R-band magnitudes we re calculated 



from t he brightness of the comparison star as in IStaude et al 
(2003). 



Simultaneously to the observations with XMM-Newton on 
March 31, 2005, 7.6 h of optical data were acquired with expo- 
sure times of 60s (see Fig. [TJ. For the period analysis, further 
data taken on March 28 and 30, and on April 20, 2005 were 

included. 

Compared to the R-band data presented by Staud e et al 



(2003), the object was almost one magnitude fainter in 2005. 
In 2003, the optical short-term variability occurred on the pe- 
riod that we would now identify with the true white-dwarf spin 
(see Sect. @), whereas in 2005 it occurred on co-Q,. 
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Fig. 3. top: R-band light curve from May 8, 2006, bottom: The 
above data with BJED<2453864.56, binned into 2 bins on the 
spin p eriod calculated with the ephemeris given by lStaude et al. 
(12003b . The error bars mark the standard deviation of the values 
entering each bin. 

There is no optical photometry simultaneous to the XMM- 
Newton observation of 2006. Nevertheless, the star was ob- 
served for 4.4 h on May 08, 2006, using the AIP 70cm- 
telescope with 60 s exposure time (see Fig. [3). 

As found for data acquired in 2003, the system was mea- 
sured to be one magnitude brighter than in 2005 and its vari- 
ability to be dominated by the white-dwarf spin period (see 
Sect.©. 

3. Analysis of the averaged X-ray spectra 

The EPIC-PN and MOS1/2 spectra for both observations are 
shown in Fig. @] Both spectra appear similar and have flux 
in all energy bins in the interval from 0.3 to llkeV. In the 
energy range below 0.6 keV, a distinct soft component is ob- 
served, which confirms the classification of MU Cam as a soft 
IP. Around 6.5 keV, contributions from Fe emission lines are 




channel energy (keV) 



Fig. 4. The EPIC X-ray spectra from 2005 (top) and 2006 (bot- 
fom)(green PN, red+black: MOS1/2) with overplotted fit and 
residuals (Parameters in Tab. |2). The insets show an enlarged 
view of the iron-K region. 



detected. The spectral analysis of our data was performed with 
Xspec 11.3.0. 



Staude et al.1 (|2006) fitted the spectrum from 2005 with 
consisting of a hot emitting gas (mekal), a 

a Gaussian represent- 



a model 

black body component (bbody), and 
ing iron-K emission lines (gaussian) absorbed by a photo- 
electric absorber (wabs), and an additional hot-gas model 
absorbed by a photoelectric absorber of different column 
density, i.e. wabs\(mekal\)+wabs2(mekal2+bbody+gaussian). 
Since the hydrogen column density, nu, of wabs i was 250 times 
higher than that of wabsi («h=2.7x10 22 cm -2 ), the model con- 
sisted of components absorbed by a hydrogen column compa- 
rable to the interstellar absorption, and an additional highly- 
absorbed mekal. 

The data from 2005 and 2006 were fitted in the range from 
0.3 to 1 1 keV with Xspec using an identical multi-component 
model for both data sets. This model consisted of an emit- 
ting hot-gas component (mekal), a black body (bbody), and 
two Gaussian emission lines representing the iron-K lines 
at 6.4 keV and 6.7 keV (gaussian), which were all absorbed 
by a partially-covering photoelectric absorber (pcfabs), and a 
simple absorber representing the interstellar medium (wabs), 
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Table 2. The main parameters of the spectral model for the 
average XMM-Newton spectra from 2005 and 2006. The num- 
bers in brackets are the errors and have to be multiplied by the 
same power of ten as to that for the parameters, 'n.h.p.' is the 
'null hypothesis probability'. 







2005 


2006 


wabs 


nH [10 22 cm- 2 ] 


4.1E-02 (0.8) 


4.7E-02 (0.6) 


pcfabs 


nH [10 22 cm" 2 ] 

V. Vll 1 CX\* L 


7.9 (0.6) 
61 (0 021 


9.5 (0.5) 
67 (0 on 


mekal 


kT [keV] 
norm 


79.9 [35] 1 
7.3E-03 (0.4) 


79.9 [35] 1 
6.5E-03 (0.3) 


bbody 


kT [keV] 
norm 


5.9E-02 (0.07) 
1.9E-04(0.9) 


5.4E-02 (0.05) 
2.2E-04 (0.9) 


gaussian 


LineE [keV] 
Sigma [keV] 
norm 


6.4 (- 2 ) 
0.1 (- 2 ) 
3.1E-05 (0.7) 


6.4 (- 2 ) 
0.07 (0.02) 
2.9E-05 (0.6) 


gaussian 


LineE [keV] 
Sigma [keV] 
norm 


6.7 (- 2 ) 
0.1 (- 2 ) 
3.1E-05 (0.7) 


6.7 (- 2 ) 
0.09 (0.03) 
2.5E-05 (0.7) 


xl (d.o.f) 




1.034659 (948) 


1.046327 (1754) 


n.h.p. 




2.2E-01 


8.7E-02 



1 : The temperature of the mekal component ran against its upper limit. 
The value in brackets is the plasma temperature which is compliant 
with both the XMM and INTEGRAL spectrum (see Fig.[5]and chapter 



3). 



: The parameters were fixed to these values. 



i.e. wabs(pcfabs(mekal + bbody + gaussian + gaussian)). The 
model fitted the data as well as the previous model, although 
it required just one partially-absorbed hot-gas component, in- 
stead of two of different homogeneous absorption. The model 
can be understood according to the standard picture of IPs: The 
bbody represents the soft X-ray emission below 0.7 keV, either 
thermal emission or reprocessed hard X-rays, the mekal com- 
ponent describes the hard component emitted by the hot post- 
shock material above the white dwarf's surface, and the gaus- 
sian emission line at 6.4keV (the Fe fluorescence line) are trac- 
ers of reprocessed X-rays. The pcfabs models absorption inside 
the system by the accretion curtain/stream. It is mainly respon- 
sible for the absorption of soft photons in the mekal compo- 
nent, which explains the relative weakness of flux in the range 
between 0.7 and 7 keV with respect to flux at the harder end of 
the spectrum. 

Values of the free-fit parameters were determined by x 2 - 
minimisation and are listed in Table|2] The fits were successful 
in both cases with a reduces xl ~ 1 ■ Adding further compo- 
nents such as another mekal or a black body, improved the fit 
only marginally. A third Gaussian was not required to model 
the Fe-line at 6.96 keV, because this line was described well 
already by the mekal model. Our fit was not influenced signifi- 
cantly by allowing the metal abundance of the mekal to vary as 



a free parameter; in this case, the mekal, in fact, hardly changed 
from its initial value. The lines of iron were resolved in data ac- 
quired in 2006 and their widths were therefore allowed to vary, 
whereas they were unresolved in the 2005 data (see insets in 
Fig. @). In the data from 2005, there is line emission at the po- 
sitions of all three iron-K lines, although an insufficient number 
of counts do not enable the lines to be resolved. We therefore 
fixed their line widths to 0. 1 keV in the fitting, assuming that 
they were as broad as in 2006. 

The parameters of the partial absorber in our fits are 
comparable to the values obtained for the other soft IPs 
observed with XMM-New ton, V405 A ur (n+ 

PQ Gem 



200. 



4), 



Evans & Hellier 

" 2 CvrFract=0A5. lEvansetafl |2006|) 



CvrFract-Q.52, 
(« H =ll.lxl0 22 cm 

UU C ol (n H =10xl0 22 cnr 2 , CvrFract =0.5 l. lde Martino et al 
1 20061). and NY Lup (n H =9.7xl0 22 cnr 2 , CvrFract=0A7, 
Haberl et al.ll2002h . 



A problem, however, with our fits is that the fitted tem- 
peratures from the mekal models are at the upper limit for the 
model available in Xspec. Reducing the temperature to lower 
values, produces fits in which the pcfabs jiH is reduced but 
its covering-fractio n increased; oervall have we a higher xl- 



Hab erl et al.1 (|2002) report that fitting a multi-temperature dis- 
tribution of mekal models with temperatures following a power 
law to the data of NY Lup, leads to an unconstrained maximum 
temperature, i.e. a similar problem. 

MU Cam was d etected to be an INTEGRAL/IBIS source. 



Barlow et al. (2006) measured a bremsstrahlung temperature of 



8.1(4.7)keV. Judging from the XMM-data, this value is clearly 
too low; we therefore analysed bot h data sets together . The 

" J2OO6J) 



was 



INTEGRAL spectrum discussed in iBarlow et al. 
kindly provided by the authors and is included in Fig. [5] Fitting 
both data sets simultaneously with Xspec produced no satis- 
fying results because of the large differences in the number of 
data points and their errors. However, by fixing the temperature 
of the meA:a/-component to a number of different values and 
fitting the remaining parameters, good fits to the XMM-data 
(xl < 1.1) and a convincing agreement with the INTEGRAL 
data points (the model was located within the error-boxes of 
the three data points) were achieved for kT = 35 + lOkeV. 

Utilizing instruments capable of observing a wider range 
of X-ray energies (e.g. RXTE), the pla sma temperatures in IP s 



can be determined more precisely (ISuleimanov et al] 120051) . 
Compared to the temperatures measured for other IPs, the value 
of fcr=35keV is unusually high. Typical values for IPs are 10- 
20 keV. Since this temperature must be regarded as an average 
of the entire temperature range, and since the reflection compo- 
nent present in the star redistributes photons to higher energies, 
the post-shock temperature remains only weakly constrained. 

The EPIC-spectra and the parameter values of the fits to 
the two data sets of 2005 and 2006 appear to be quite similar, 
which is surprising when differences, particularly in the soft 
X-ray light curves are considered. 

The spectra from RGS 1 and RGS2 of both epochs were 
grouped to contain at least 20 counts. The fit to the integrated 
EPIC spectra reproduces the data well. No significant lines 
were detected in the wavelength range considered, and we 
therefore disregard the RGS data in further analysis. 
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MU Cam 



r: 




0.001 0.01 0.1 1.0 10.0 100.0 

Energy [keV] 

Fig. 5. The spectral energy distribution of MU Cam. Shown 
are the data synchronously obtained in 2005 (rhombs and X- 
ray model spectrum in black), the data from 2006 (XMM- 
Newton: circle and X-ray model spectrum in red, /?-band 
flux from the 70cm telescope: triangle), and, additionally, 
2MASS J, H, and K fluxes {squares) plus the opti cal spec 



tra fro m lArauio-Betancor et al.l d2003l) {lower) and IWei et al 



19991) {upper). Additionall y the three data p oints from the 



INTEGRAL/IBIS spectrum Barlow et alj 120061) are plotted 
with their corresponding error bars. The X-ray models are plot- 
ted without the interstellar absorption component. To illustrate 
the shape of the black body component in the X-rays, the ex- 
trapolation to lower energies outside the fitted range is plotted 
with dots. 



The overall spectral energy distribution of the system is 
shown in Fig. [5] where the available measurements are plotted 
with colours that differentiate between data obtained simulta- 
neously (black and red), acquired in an optically-bright state as 
assessed from the R-magnitude (blue), and of unknown state 
(green). 

During an epoch of high optical brightness, the X-ray flux 
is lower than when the system is optically faint. It is evident 
that the black body component, which can be seen in the X- 
rays, contributes only negligibly to the optical/UV flux. The 
shape of the continuum in the optical/UV wavelength range 
resembles that of black body with a much lower temperature 
and can be attributed to the accretion disk. 



4. Variability analysis 

The light curves derived from the observations during the two 
epochs are shown in Figs. Q] and 12 the corresponding peri- 
odograms in Fig. [6] The analysis is supported by additional 
optical photometric data that are listed in Table Q] In spite of 
differences between data for different epochs, the periodograms 
provide a clear measurement of the basic periods of the system. 



4.1. The basic frequencies 

Arauio-Betancor et all d2003h and Istaude et all d2003l) gave a 



ods from optical photometric data. Op tical spectroscopy was 
used by lArauio-Betancor et al.l q2003) to identify the orbital 
frequency, SI. Using the optical data alone, it was impossible 
to identify unequivocally the white dwarf spin frequency, co. 

Figure [6] shows the results of a p eriod analysis using the 
meth od of analysis-of-variance (AOV, ISchwarzenberg-Czernv 
1989) in the different energy bands, for data from the two 
XMM campaigns. We indicate the orbital frequency Q, the ten- 
tative spin frequency co, the side-bands of both these frequen- 
cies, and the (sub)harmonics. AOV tends to create artifacts in 
the form of sub-harmonics; this is in contrast to the Fourier- 
transform, which creates harmonics. The main peaks in the X- 
ray bands HX and MX are in both datasets at the assumed co. In 
the SX band, the pattern is more complex in 2005 than in 2006, 
showing variability not only at co but also at Q and 2co-Q. The 
OM-data show variability in the so-called beat-frequency co-Q, 
at least in data from 2005. In the optical R-band, the peak is 
located at either co (2003 and 2006) or co-Q (2005). 

A consistent explanation for the presence of all peaks in 
the periodograms can be given when the tentative value of co is 
assumed to measure the white dwarf spin. 

The X-ray periodograms from 2006, which show variabil- 
ity power only at Q, co, and 2co, could be produced by accre- 
tion from a truncated disk to two opposite acc retion regions 
via shock-heated curtains dWvnn & Kinglll992l) ; reprocession 
on the white dwarf surface would then produce the soft X-rays 
{SX). In 2005, the harder X-rays {HX and MX) were created by 
accretion along curtains, indicated by variability in mainly Q 
and co. The soft X-rays (SX) have variability power also in the 
frequency 2co-Q and some signal at co-Q.. This can be explained 
by their origin in reprocession on the white dwarf as well as by 
accretion directly from the accretion stream via magnetic field 
lines, wh ich would produce m odulations in superpositions of 
Q and co dWvnn & Kinglll992h . 

In 2005, the UV data shows variation mainly at co-Q,, the 
periodogram has weak power only at Q and co. The 7?-band 
data has a similar variability pattern with stronger variation in 
Q. Changes in the dominating frequency from co to co-Q, es- 
pecially at the lower energies, can be described by changing 
accretion modes within the current models of IPs. 

We could imagine assuming the tentative co to be the true 
co-Q. This would be consistent with diskless accretion along 
magnetic field lines. Certain problems arise in this approach, 
for example the need to explain the absence of any signature 
of variation with frequency co. Since the visibility of the accre- 
tion regions is modulated by this frequency, it is expected to 
be observed in hard X-rays for which absorption is probably 
incomplete. In the /?-band and UV data from 2005, one would 
have to explain the strong signal at co-2Q (the tentative co-Q), a 
small contribution to co-Q, and no trace of co. 

Warner (1986) explained an co-2Q varibility to be an or- 



bital modulation of the co-Q component. It is therefore unlikely 
that the latter has a far smaller amplitude. We would also have 
to replace the conv enient explanation of 2co-Q proposed by 



Wynn & King! (119921) . by one for 2co-2Q; furthermore, for the 



prelimininary identification of the system orbital and spin peri- 



2006 data, we would have to consider 2co-Q without a signature 
for the 2co periodicity. 
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7? -band photometry - for comparison, the AOV analysis of the 7?-band data from 2003 ( Sta ude et al.ll2003l) is added; right: The 
data from 2006 (Fig.[2j in different energy ranges and instruments. The orbital (£2) and spin (a)) frequencies, their sidebands, and 
(sub-)harmonics are marked with dotted lines. 



To summarize, it is unlikely that this second interpretation 
reproduces our observations of MU Cam. We th erefore claim 



that w e have found the spin period predicted by Staude et al 
(2003) to be the true rotation of the white dwarf. 



Aside from their common properties, the data from the two 
epochs show a number of differences. 

In the light curves from 2005 (Fig. [T}, the soft X-rays (SX) 
and the R-band data show a strong dependence of variability 
amplitude on the orbital phase. In the soft X-ray light curve, 
two peaks per spin-cycle of various heights are directly visible. 
They are probably the origin of the frequency 2a>-£l. 

The most striking differences between the light curves ob- 
tained on April 06, 2006 (Fig.©, and those from 2005 (Fig.Q), 
that there is a reduction of the SX countrate and the ampli- 
tude of the short-term variation does not depend obviously on 
the orbital period. The orbital period Q is detected in the peri- 
odograms of all the X-ray bands, although there is no obvious 
variation of this timescale in the light curves (see Fig.lTOt. 

The R -band data - in agreement with the XMM-Newton 
data, although these were acquired one month later - demon- 
strate variability mainly with u>. Due to the short duration of the 
observation - slightly less than one orbit - no secure statement 
about the variability on the orbital period can be made. 



4.2. Towards a refined spin ephemeris 

Using the data acquired in 2005 and 2006, we attempt to im- 
prove the accuracy of the spin ephemeris. The optical data of 
long-scale base line can be used to refine the period, while the 
hard X-rays relate the ephemeris zero-point to the rotation of 
the white dwarf. 

From the spin-phase folded light curve of the R-band data 
from 2006 (Fig. [3J, a new time of spin maximum for a long- 
term ephemeris can be derived. The phase of maximum, 0.075, 
corresponds to 7 , max sp i n (BJED)= 2453864.5668(3). The error 
in the epoch is assumed to be half the bin width, 0.025 spin 
phases. The cy cle-count is lost du e to the error in the spin 
ephemeris from lStaude et al.l (12003I) . 



Kim et alj d2005l) proposed a refined spin ephemeris, 



p jpi „=0.01374116815(17)d, based on a 3-period fit. When 
folding our R-band data from May 8, 2006 with their 
ephemeris, the spin maximum, however occurs at phase 
0.725(25), which is clearly outside the error ranges of their 
ephemeris zero point at the time of observation (0.01 phases). 
The reason for this discrepancy may be that they used data ob- 
tained in different accretion states: this is obvious from the sig- 
nificant change in brightness between their measurements in 
2004 (m R ~ 14.5) and 2005 (m R ~ 15.5). We detected fun- 
damental differences in the variability patterns under similar 
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conditions, i.e. a change from variation with frequency a> to 
frequency oj-D.. 

We calculated an X-ray spin-ephemeris from phase- 
binned light curves of the harder photons (MHX, 700 - 
12000eV, Fig. |7). Times of maximum flux were de- 
temined from a spline fit to the binned data to be 
7max,PN haill (BJED) = 2453 461.4475(8) (phase 0.68 +0.02) and 
7 , max,PN hanl (BJED) = 2453 832.0615(8) (phase 0.57+0.02). The 
assumed uncertainty in the maximum-determination is one bin- 
width. The flux maxima were used because they represent the 
time of optimal visibility for the principal accretion region of 
the white dwarf. Depth and position of the minimum may be 
strongly dependent on the size of the accretion region(s), be- 
cause their extension determines the length and completeness 
of their occultation behind the white-dwarf limb. The apparent 
optima visibility of the main accretion region is only weakly 
dependent on the size. 

(12003b provides a 



The spin-ephemeris of IStaude et al 



cycle-count of 26971, which corresponds to a period of 
P« pWffl r=0.01374120(6)d. 

This value matches the previously-published results well. 
The measurement of period confirms the previous estimate for 
an independent data set, even though: (i) the period does not 
appear to be more accurate than that derived from o ptical pho 



tometr y in 2003; and (ii) the spin period measured bv lKim et al 
(2005), is within the error range of our results. 



The timing information derived from optical photometry 
cannot be combined with that of the X-ray data, because the 
two types of radiation originate in completely different parts 
of the binary system, even if their measured variabilities have 
similar periods: a time-lag between the pulses is therefore ex- 
pected. 

4.3. From periodicity to structure 

In 2005, the spin-folded light curve in band MHX is sawtooth- 
shaped (Fig. |7), with a slow rise to maximum in ~0.7 phase 
units, while the decline lasts only 0.3 phase units. The light 
curves have a pulsed fraction of 50 per cent and show a similar 
behaviour in the energy sub-bands MX and HX. 

In 2006, the MHX light curve looks similar to that from 
2005, but its minimum is more pronounced with a pulsed frac- 
tion of 75 per cent compared to 50 per cent in 2005. 

To obtain more information about the origin of the spectral 
components, a phase-resolved spectral analysis was performed 
with respect to different periodicities (a>, £2, w-Q, 2w-Q). The 
data were binned into 10 phase bins for each period. A num- 
ber of parameters from the model fit to the global spectra were 
fixed to the values previously obtained: wabsjiH (since it is in- 
terpreted to be of interstellar origin), pcfabsjiH, mekalJcT (be- 
cause this was, in any case, not constrained), bbodyJcT (which 
when allowed to vary, changed very little around its mean 
value), and gaussianf l,2)jsigma. 

The error bars in the binned light curve plots indicate the 
standard deviation of values that enter the bin. We calculate the 
errors in our phase-binned dataset by measuring the Poisson 
noise of photons in a given bin. This is the correct approach 
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Fig. 7. The background-subtracted PN count rate in the en- 
ergy band MHX (700 to 12000eV) (top 2005/03/31, bottom 
2006/04/06) binned to 30s, plotted over spin-phase (upper). 
The same data, binned into 50 phase bins with an overplotted 
spline fit (lower). 



when a system is observed repeatedly in a similar state, which 
is clearly the case if data show variability only in the period for 
which the data are binned. 

When binning data that are varied by multiple periods, the 
measurements from different cycles, entering the same bin, are 
not in general obtained when the system is in a similar state. If 
Poissonan noise was assumed to represent the error, this would 
not take account of the intrinsic variability by the superposition 
with the other periods, and would underestimate the errors. 

In our data for MU Cam there is always variability on mul- 
tiple periods - in the simplest case for just to and Q. We could 
attempt to remove the variation that occurs on one period from 
the data, but this will certainly not produce the desired result 
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Fig. 8. The spin-folded PN light curves from for 700-12000eV (MHX) and 200-700eV (ST). For the calculation of the spin- 
phase, we u s e as z ero points the epochs of maximum flux in the MHX band determined in Sect. 14.21 and the period from 
Staude et al.l ( 12003b . The error bars in the light curves mark the standard deviation of the values having entered the bin. Below 
are plotted the spin-phase dependent parameter values for the fit to phase-selected MOS+PN spectra. The parameters not shown 
here were kept constant at the values determined from the fit to the average spectrum, left: March 31, 2005, right. April 06, 2006 



in the presence of strong side-bands, as in the SX-band in the 
2005 data. 

To maintain the size of errors in all plots of binned light 
curves, we consistently use the standard deviation of values as 
they enter each bin throughout the paper. These values repre- 
sent upper limits to the true errors. 



4.3.1 . The white dwarf rotation - to 

The parameter values of the fits to the spin-phase selected spec- 
tra are shown in Fig. [8] with phase-folded light curves in the 
bands SX and MHX . The similar shape of the soft phase-folded 
X-ray light curves is remarkable, despite the obvious differ- 
ences in the time domain (compare Figs. [TJ and It appears 
that the strong variations in the 2005 data are completely re- 
moved when the data are folded and averaged using the spin 
period. 



4.3.2. The stars' orientation matters - oj-Q. 

Variability at the beat-period w-Q was detected only in the data 
from 2005. It was the dominating period in the OM and the R- 
band light curves, it was clearly present in the soft X-ray band 
(SX), whereas it was almost undetected in the hard X-ray bands. 
The beat-phase folded X-ray light curves and phase-resolved 
spectra do not show significant variation and are therefore not 
shown. 



4.3.3. A signature of two-pole accretion - laj-Q. 

While the X-ray data from 2006 varies predominantly at the 
spin period, the soft X-ray component in 2005 shows varia- 
tion also with additional frequencies, most prominently 2w-Q. 
This is an indication that there are two accreting poles, whose 
accretion rates are dependent on the orientatio n of the white 
dwar f within the binary frame, i.e. of both stars ( Wvnn & K ing 
1992b . To decide whether there is any difference between the 
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emission of the two accretion spots, the soft and hard X-ray 
photons were phase-folded on the period corresponding to (±>- 
Q/2. Spectra were extracted for ten phase-bins with respect to 
this period (O w _ n/2 = (BJED - 2452682.4181)/0.01423886) 
and fitted. The light curves and the fitted parameters are shown 
in Fig. i 

Separated by 0.5 phases, there are two humps that have sim- 
ilar a soft X-ray light curve, which explains the strong signal 
of frequency 2a>-Q.. The parameters of the fit indicate that the 
source of this variability are two distinct poles: the mekal com- 
ponent has one peak only above an almost constant level, and 
the integrated black body normalization of the humps clearly 
differ. 
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Fig. 9. Light curves in the hard and soft X-ray bands from 2005 
folded on the period corresponding to co-Q/2 (error bars in the 
light curves are the standard deviation of the values having en- 
tered the bin), and the free parameters of the fits to the phase- 
resolved spectra. 



4.3.4. The orbital period - Q. 

For the data discussed here, a significant and likely periodic de- 
pendency of the radiation on the orbital motion is detected only 
for 2005 (see Fig.|6]l. Compared to the power in other detected 
frequencies, the orbital signal is strongest in the soft X-rays, 
and detected clearly in the other X-ray bands and the optical. 
Phase-folded light curves and fitted parameter values for the X- 
ray spectra from 2005, are shown in the left panel of Fig. [10] 
There is almost no variation in the UV light curve (which is 
therefore not shown), but the other light curves display clearly 
a single hump. The maxima and minima in the curves do not 
occur at the same orbital phases - the R-band is 0.1 phases 
ahead. The probable origin of the variability CI is the changing 
visibility of emission sites during the rotation of the binary. 

In contrast to the spectral fits applied to data averaged us- 
ing the other frequencies (see previous sections), the parameter 
pcfabsjiH was allowed to vary for studying the Q-dependence 
of the data from 2005. By fitting the binned data, we searched 
for variations in the hydrogen column density, caused by the 
changing orientation of the magnetically-guided stream with 
respect to the line-of-sight during one revolution. The number 
of observed orbital cycles, however is rather low for the binned 
data to be regarded as representative of the mean system; since 
fewer than two consecutive orbits were observed, the data can 
be considered only as a snap-shot of the system at a certain 
time. 

The data from 2006, phase-folded on Q, are not shown, be- 
cause only marginal variability is observed in the X-ray light 
curves and OM data; the latter data set and optical data are 
affected by insufficient coverage of the orbital cycle. We there- 
fore do not attempt any phase-dependent analysis. 

5. Discussion 

5.1. Pure disk accretion observed in 2006 

In all of the X-ray periodograms, the dominant frequency is the 
white-dwarf spin <d. The frequency 2u> is also present, but far 
weaker. In the periodograms, there is signal at the orbital fre- 
quency Q. in all wavelength ranges, although it is not prominent 
in X-rays and not completely resolved in opti cal light. 

Following the analysis of IWvnn & King ( 1992b . the natu- 
ral explanation for the detected variability pattern is accretion 
via an accretion disk and accretion curtains, which connect the 
disk to the magnetic poles. The weak detection of signal at 
frequency 2co, in addition to that at to, implies that there are 
two accreting spots, almost opposite to each other Their prop- 
erties or visibility must differ to explain the far stronger signal 
at u>. This interpretation is supported by Fig. [8] (right panel), 
where the soft X-ray light curve shows (when compared to the 
hard profile) a second maximum shortly before phase 0.5, and 
the normalisation of the black body component of the phase- 
resolved spectral fits has a secondary maximum at this posi- 
tion. The normalisation of the me&aZ-component shows only 
one maximum - located at the phase of maximum brightness - 
with a slow rise and a rapid decay. The sawtooth-like shape in- 
dicates its origin in an arc-shaped region. The covering-fraction 
of the absorption component in the spectral fits shows a smooth 
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Fig. 10. The orbital period Q: Light curves from 2005 in the 
hard and soft X-ray bands, the fitted spectral parameters, and 
the phase-folded /?-band data. The error bars are the standard 
deviation of the values having entered the bin. 



modulation on the spin period with a clear maximum at the 
phase of minimum X-ray flux. 

The resemblance of the spin-folded X-ra y light curves in 



Evans & Hellier 



the tw o energy ranges to the ones shown by 
d2004 for V405 Aur is striking. As observed there, the soft X- 
rays show a double-hump, whereas the hard ones form a saw- 
tooth shape. A similar accretion geometry, therefore, appears to 
exist in MU Cam, where the white-dwarf dipole axis is located 
close to the orbital plane instead of being aligned with the spin 
axis. The hard and soft X-ray light curves and the parameters of 
the fits to the spin-resolved X-ray spectra (Fig. [8] right panel) 
can be explained in a model where the lower accreting pole has 
optimal visibility at spin phase 0.45, and the upper at phase 0. 
Between phase and phase 0.3, the covering fraction of the 
partial absorber increases. This absorbing matter is the accre- 
tion curtain to the upper pole, crossing the line-of-sight to the 
X-ray emitting region. The phase-shift between the maximum 
absorption and the optimal visibility of the accretion region can 
be explained by the bulk of matter being accreted along trailing 



field-lines, as proposed by iBeardmore et al.l dl998l) to explain 
absorption features in observations of FO Aqr. At phase 0.25, 
the shocked region above the upper pole, which emits hard X- 
rays, disappears behind the white dwarf, in a similar way to the 
soft X-ray spot. An indication that soft X-ray emission traces 
reprocession sites of hard X-rays is the modulation of soft X- 
rays by the white-dwarf spin. 



In contrast to the optical data discussed in lEvans & H ellier 



(120041) . our R-band data from 2003 and 2006 do not show a 
double-peaked spi n pulse. The spin- folded optical data from 
2003 (see Fig. 4 in lStaude et a02003l) have a zigzag-like shape 
of amplitude ~0.2mag. The spin-phase folded R-band light 
curve from 2006 (Fig. [3]) shows a sinusoidal variation of ampli- 
tude 0. 1 mag. Using the hard X-ray determined spin-ephemeris 
(Sect. |4~2l , the flux maximum for the 2006 data occurs at phase 
0.55, i.e. half a spin cycle later than the X-ray maximum. We 
propose that the pulsed fraction (10 per cent) of the optical flux 
to be originates in the irradiated parts of the accretion curtains. 
Their irradiated inner regions are most easily visible when the 
upper pole points away from the observer. The magnetic axis 
of the white dwarf must be tilted, with respect to the rotational 
axis, by approximately the system inclination angle to be able 
to explain the maximum X-ray flux that is half a spin-cycle 
offset from the optical maximum. When pointing towards the 
observer at phase zero, this prevents obscuration of the upper 
accretion region by its accretion curtain. 

Unfortunately, the OM light curve shows a clear variation 
pattern only towards the end of the observation and is affected 
by data gaps; it is therefore impossible to cross-correlate the 
OM light curve with the X-ray and optical light curves to mea- 
sure the relative phasing, i.e. the position of the reprocessing 
sites. 



5.2. Accretion via disk and stream observed in 2005 

The hard X-ray component is assumed to be originate in a 
shock inside the accretion column above the magnetic pole. 
Since there is almost no variability in the count-rate of the hard 
photons on the beat frequency co-Q., the accretion rate does not 
appear to be dependent on the white-dwarf orientation within 
the binary frame. The accretion producing the emission of hard 
photons probably happens via an accretion disk, as occured in 
2006. 

The number distribution of the soft X-ray spectrum (SX, 
200 - 700 eV) is reproduced well by a black body, and appear to 
originate on the white dwarf surface (because of the strong sig- 
nal at co). The presence of co-Q. shows that the orientation of the 
white dwarf in the binary frame has influence on the emission. 
Stream-fed accretion, penetrating and heating the white-dwarf 
surface, is probably therefore the origin of (parts of) the soft 
X-ray component, because pure reprocession of hard X-rays 
on the white dwarf would produce a signal at co only. The fre- 
quency 2<y-Q implies that there is a second accreting pole that 
is opposite to the principle pole, since it denotes the detection 
of a similar accretion geometry after just a little more than half 
a spin -cycle, in the case of stream-fed accretion (iWvnn & K ing 
19921) . The two peaks per spin cycle can also be directly ob- 
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served in the light curve (Fig.[T). Figure [9] shows that only the 
soft X-ray flux (i.e. the flux in the black-body component), not 
the hard X-ray flux, is modulated with respect to the orientation 
of the two stars. 

As the dominant frequency in the optical and UV light 
in 2005, the beat-frequency oj-D. indicates that the emission 
or visibility of this radiation component was strongly depen- 
dent on the orientation of the white-dwarf magnetic field in 
the binary reference frame. The accretion disk may have been 
smaller in 2005 than in 2006, and the inner rim therefore lo- 
cated at a greater distance from the white dwarf. This would 
allow a disk-overflowing part of the ballistic stream, which 
would otherwise be negligible, to dominate accretion via the 
disk. Only the prolonged ballistic stream would have a suffi- 
ciently high density and incident flux to re-emit significantly, 
producing a co-Q. periodicity. This view is supported by the 
lower brightness of the system in the optical/UV in 2005. It 
is observed that the amplitude of the short-term variations did 
not vary significantly, in contrast to the size of the constant 
component. 

The orbital variation in the X-ray emission (Fig. ITOb is 
probably caused by the different visibility of the accretion sites 
in the binaries frame of reference. Our fitted parameter s argue 
against a model, pushed forward by lHellier et al.l (1199 lh to ex- 
plain orbital modulation of the X-rays in AO Psc, where the 
minimum at orbital phase 0.3 is explained by increased absorp- 
tion, e.g. by matter lifted out of the orbital plane. Unfortunately, 
the published orbital ephemerides either have to o large errors 



the published orbital ephemeriaes either have too large errors 
( Araujo-Betancor et alJ 2003[ Staude et al.l 2003) or a remain- 
ing cycle-count ambiguity (IKim et al.ll2005l) to allow to relate 
the relative phasing of the other detected periods to the proba- 
ble epoch o f inferior conjunction of the s econdary star, as de- 
termined by Araujo-Bet ancor et al. (2003) from phase-resolved 
spectroscopy. 

To explain the behaviour of the system during the XMM- 
Newton observation in 2005, we therefore propose that the 
white dwarf was accreting via an accretion disk and a disk- 
overflowing stream. The hard X-ray emission, which is partic- 
ularly sensitive to variations on the white-dwarf spin period, 
are produced by shocks above the footlines of the accretion 
curtains being fed by the accretion disk. The main (upper) ac- 
creting pole is best visible around spin phase 0, when it points 
to the observer. The occurrence of the frequency 2co-Q. in the 
soft X-ray band is an indicator of two opposite, similarly bright, 
accretion spots, accreting via an accretion stream. This is evi- 
dence of the white dwarf's dipole axis being strongly inclined 
from the rotational axis, which produces similar visibility and 
accretion efficiency in both regions just offset by 0.5 phases in 
spin and beat period, respectively. 

In the spectral model, we included two Gaussians to rep- 
resent iron lines at 6.4 and 6.7 keV. Although they have no 
significant influence on the overall they may support the 
interpretation of the spin-dependent behaviour of the system. 
The normalisation of the Fe -line at 6.7 keV closely follows the 
shape of the normalisation of the black body component with 
two distinct peaks. The fluorescence line at 6.4 keV appears 
to have only one peak and an eclipse-like feature at spin-phase 
0.7, which resembles the mekal component. In general, this line 



is interpreted to be caused by the reflection of the hard X-ray 
component by cold material, e.g. by the irradiated surface of 
the white dwarf close to the accretion region. The coupling of 
the normalisations of the mekal to the 6.4 keV line supports this 
view in the case of MU Cam. 

In the optical spectra published for MU Cam dWei et al 



19991 I Araujo-Betancor et al.l2003l) . there is no evience of spec- 
tral features that could be produced on the secondary star. It is 
therefore impossible to estimate the contribution of the late- 
type star to the overall flux and the distance of the system. 
Using orbital period-spectral type relations for cataclysmic 
variab les determined either observationallv ([Smith & Dhillonl 
19981) or theoretically (IBeuermann et alJI 19981) . the companion 
star in MU Cam is expected to be of early M type (i.e. ~M1). 



6. Summary and conclusions 

We have presented the analysis of comprehensive data sets ob- 
tained simultaneously at optical, ultraviolet, soft and hard X- 
rays in 2005 and 2006 with the XMM-Newton observatory and 
the AIP 70cm telescope. Our main results and conclusions may 
be summarized as following: 

- Th e period ten t atively derived as the white dwarf spin 
by IStaude et alJ d2003l) is the dominant signal in the hard 
X-ray light curves. We claim that this provides a mea- 
sure of the white-dwarf rotation. By linear interpolation 
of the timings of the mean hard X-ray pulse maxima in 
the two XMM-Newton observations we measure a spin 
ephemeris of B7£D(max, PN hard ) = 2453 461.4475(8) + 
Ex 0.01374120(6) 

- A distinct soft component in the X-ray emission, which was 
already visible in the ROSAT spectrum, was unambigu- 
ously detected and could be fitted well with a black body 
of temperature 59 eV and 54 e V in 2005 and 2006, respec- 
tively. MU Cam could be reliably identified as a soft IP. 

- MU Cam shows changes in its accretion state between disk- 
dominated and disk with an additional stream component 
during the observing interval from 2003 to 2006. These 
are displayed by a change in the mean optical brightness 
of about 1 mag, and by fundamental changes in the power 
spectra throughout the complete observed energy range. 
Such changes were also reported for other IPs (e.g. TX 



Col: Norto 



Evans et al 



et al.Ml997L FO Aqr: iBeardmore et all 1998 



2004), and therefore appear to be a prevalent 



feature of intermediate polars. 

- The occurrence of the lu frequency in the X-ray data from 
2006 and the strength of 2a>-Q. in 2005 suggest two accre- 
tion regions at opposite positions with similar properties 
and visibility, at least in certain energy ranges. This may be 
achieved if the dipole axis of the magnetic field of the white 
dwarf was strongly inclined away from the rotation axis to 
the orbital plane. 

- The strong modulation in all wavelength ranges implies a 
rather high inclination, whereas the absence of an X-ray 
(white dwarf) eclipse limits it to < 75°. It is probably even 
lower, because no eclipse is detected in the optical, which 
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would be a sign of a partial obscuration of the accretion 
disk. 

- To improve our understanding of the properties of the 
system, future observations should include more optical 
photometry to refine the spin and orbital ephemerides. 
For studying the accretion geometry, time-resolved opti- 
cal high-resolution spectroscopy will be crucial, and it may 
help to determine an accurate orbital ephemeris. To de- 
termine the distance, near-infrared spectroscopy should be 
performed to be able to determine the spectral type of the 
secondary star and its contribution to the spectral energy 
distribution. 
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